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ABSTRACT 


Simple formulas for calculating reflection and transmission coeffi- 
cients for converging-diverging ducts are derived; they extend the method of 
Cho and Ingard to arbitrary, slowly varying ducts. These formulas involve 
two parameters: The first is a function of duct snape and the second is the 

ratio of the duct radius downstream of the throat to that upstream of the 
throat. An extension of the method to include mean flow is made for sym- 
metric ducts. 


There are three possible behaviors of acoustic waves propagating in 
converging-diverging ducts: 

(1) They may propagate freely through the duct. 

(2) They may propagate in the section upstream of the contraction and 

be cut off downstream, hence suffering total reflection at the contraction. 

(3) They may propagate upstream and downstream of the contraction, but 

be cut off within the contraction, hence suffering partial reflection while 

oeing partially transmitted. It is the third case that is of concern in 
this report. 

This case has been treated by Cho and Ingard (ref. 1) for waves in a 
circular-cosh duct. The present report extends their evaluations of reflec- 
tion and transmission coefficients to nonuniform ducts that are not of 
circular-cosh shape. Our approach is to use the propagation equation of the 
circular-cosh duct as a comparison equation with one free parameter. An 
"eigenvalue" relation is derived for that parameter, and it, in turn, deter- 
mines the reflection and transmission coefficients. 

Our ultimate objective is to provide formulas for computing reflection 
and transmission coefficients for an arbitrary (slowly varying) duct. Such 
formulas are provided by equations (6) and (8). These equations require the 
duct shape as input, and after simple approximate or numerical integrations 
they yield the required coefficients. In an appendix the present method is 
extended to include a mean flow through the duct, although this can be done 
only if the duct is symmetric. The mean flow increases transmission. 

Discussions with Dr. Y.-C. Cho of Lewis have been of great help in 
writing this report. 


to 


INTRODUCTION 


SYMBOLS 


A 

b 

B 

c 


amplitude function, defined by eq. (3) 
duct radius 
defined by eq. (2a) 
sound speed 



k wave number 

t defined below eq. (13) 

M Mach number 

R reflection coefficient 

s axial coordinate of circular-cosh comparison duct 

T transmission coefficient 

U mean flow 

X axial coordinate of duct 

a duct-mode eigenvalue 

3 free parameter determining throat radius of circular-cosh duct 

Y defined by eq. (13) 

6 phase shift, defined by eq. (10) 

e maximum wal 1 slope 

ic,x reduced wave numbers, defined by eq. (5a) 
o defined below eq. (11) 

T downstream radius of duct 

p acoustic wave function for duct 

X acoustic wave function for circular-cosh comparison duct 

Subscripts: 

upstream variable, x -► -» 

0 variables at throat, x = 0 

+ downstream variables, x -► » 

The propagation equation for the n^*l-radial, m^^circumferential 
mode in a circular duct with walls of small maximum slope is 

d2^/dx2 + [k2 - a^/b^(ex)]() = 0 (1) 

(eq. (18) of ref. 1), where k is the wave number, b is the duct radius, 
o is the n^*^ zero of J|ti(o), x is the direction of propagation 
(axial direction), and e « 1 is proportional to the maximum wall slope 
(fig. 1). Equation (1) applies in the absence of mean flow. In the 
appendix this analysis (for symmetric ducts) is extended to include a one- 
dimensional gas dynamics flow through the duct. 

For the circular-cosh duct, equation (1) becomes 

d2x/ds2 + [k2 - o2/b2(gs)]x = 0 (2) 


where 


B"^(cs) = ^ (1 + t"^)( 1 + B^sech^(es))+ ^ (t"^ - l)tanh(es) (2a) 

For clarity, dependent and independent variables are denoted by x(s) and 
s in equation (2), although equation (2) is just a special case of equation 
(1). Here the duct radius is unity upstream (s ♦ -») and t down- 
stream (s -*■ *) . The parameter 3 determines the shape of the circular- 
cosh comparison duct within its contraction section. (3 is the free 
parameter alluded to in the Introduction.) 

Equation (2) can be solved exactly by putting it into the form of a 
hypergeometric equation [ref. 1 and ref. 2, section 12. 3J. However, for pre- 
sent purposes that solution is not required explicitly. It suffices to de- 
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note it by x(s). Equation (2) is used as a comparison equation (ref. 3) 
for equation (1), and an approximate solution to the latter is sought in the 
form 


t> = A(ex)x(s(cx)/e) (3) 

Substituting equation (3) into equation (1), using equation (2) to elimi- 
nate x"» snd equating to zero the coefficients of like powers of e gives 

(s')2[k2 - a2/B2(s;s)] = [k2 - a2/b2(ex)] (4) 

at 0(e0). Here primes indicate differentiation with respect to argu- 
ment. The equation at 0(e) determines A(ex), but it suffices here to note 
that A(±») = 1. Thus the x-dependence of A does not affect reflection 
and transmission coefficients and so need not be evaluated. 

Because we are dealing with case (3) of the first paragraph of this 
report, there are points and xb (resp. $a and $b) at which 
the right side (resp left side) of equation (4) vanishes. If equation (4) 
is integrated as 


where 


tc(s;8)ds = /•ex x(x)dx 


/ = [k^ - o^/B^(s;e)] 


= [k^ - a^/b^(x)] 


(5) 


(5a) 


an implicit equation for s(ex;s) is obtained with the parametric dependence 
on 8 noted. Now, for x^ < x < X5 and s^ < s < sp both sides 
of equation (5) are imaginary because r 2 and ^2 are negative when 
X is in the interval Xa to xp and s is in the interval Sa to 
sp. For X > xp and s > sp both sides of equation (5) are complex 
and can only be equal if they have the same imaginary part. Equating these 
imaginary parts gives the eigenvalue relation 


JT Sj^(8) 


ic(s; 8)ds 



x( x)dx 


( 6 ) 


The free parameter must be adjusted to make the integrals in equation (6) 
equal. The functions Sa(S) and sp(8) in equation (6) are the zeros 
of . Using equation (2a) and equating equation (5a) to zero gives 


tanh(s^^b) 


(1 - t") 

26'^(1 + - f ) 


^ (23*^) 




+ 48 


1 + 8 - 


2 k" 


a (1 




3 


( 7 ) 


Equation (6) with SglB) and 5^(6) given by equation (7) provides 
the required equation for s. The integral on the right side can be eval- 
uated for any particular duct of known cross-sectional area. The expression 
on the left side then can be inverted (numerically, in general) to find e. 

REFLECTION AND TRANSMISSION COEFFICIENTS 


The reflection and transmission coefficients for the circular-cosh duct 
are given in references 1 and 2 as 

COSh[(k^ - k_)ir/€] + COSh(oir/e) 

^ “ COShL(k+ + k_)ir/ej + COSh(oir/e) 


T = 


2 sinh(k^ir/e)sinh(k_ir/c) 
COShL(k^ + k_)ir/ej + COSh(ait/e) 


( 8 ) 


where 


kj = - .2,-2 


k2 = k2 - .2 


(8a) 


.2 = 2(1 ♦ ,- 2 ) _ ,2 ^ , ^,^2 , 


In the second form of a^, kQ is defined by 


k ^ = ^ (1 + ? (^2 + k ^) 

Clearly, k^, k^, and k^ are i evaluated at and 0. 

The value of b from equation (6) substituted into equations (8) and (8a) 
gives R and T. 

The phase shifts that occur upon reflection and transmission of the 
wave can also be determined. Letting x ± ® in equation (5), we find 


ex + s, - ex, + 

a a 


S - < 


CX + Sjj - exjj + 




a a 


ds = cX - e6_; X ♦ -® 

( 9 ) 

ds = ex - e«^; x ♦ • 
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The integrals appearing in the definitions of 6+ and 6_ have been 
put into forms that converge at the upper limits. The upstream and down- 
stream limits of the wave function (eq. (3)) are (ref. 2) 


ik s/e 


|/nr 




-ik_{S/e-6n} 

e ; X 


( 10 ) 


V. 


e 


ik^ {s/e+6-p} 


; X -► a. 


The first limit corresponds to the sum of the incident and reflected waves; 
the reflected waves are phase shifted with respect to the incident wave by 
k_6(^. Similarly, the downstream limit corresponds to a transmitted wave 
with a phase shift of k+6j. Substituting equation (9) into equation 
(10), we find that in a non-circular-cosh duct the reflected wave suffers 
an additional phase shift of 2k_6_, and the transmitted wave suffers a 
Shift of k+(6+ + 6_) - both measured with respect to the incident wave. 


SYMMETRIC DUCTS 


For the symmetric duct x = 1. With that value equation (6) becomes 


k 


2 

0 


IT 



U - / ) 'iy 

[l ♦ (kj/k^)/] 




( 11 ) 


where o = 



IX I dx 


is determined by the actual duct shape and 


hence is assumed to be known. The left side of equation (11) is obtained 
by substituting equation (2a), with x = 1, into equation (6). 

Now, if k_ (which equals k+ for symmetric ducts) is 0(1) and 
ko is also 0(1), then R « 1 and T « 0 fbr small e; that is, the 
wave is totally reflected. However, if k_ = 0(1) and kg = 0, 

T = 1/2 and R » 1/2 (ref. 1). In general, with k_ = 0(1), both 
reflection and transmission occur only when kQ/k_ « 1. Present inter- 
est is in this nontrivial case. 

These considerations suggest that equation (11) be expanded in powers 
of kQ/k_: 


e 


k2 

"o 



r(n + 2) 


(11a) 
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As asymptotic expansion for ko(«) results from inversion of equation (lia): 


4 -- 

0 ir 


k « + 


4 J- 



( 12 ) 


Because of equation (12) we have written o as a function of 
below (eq. (8)): The e dependence of R and T appears implicity 

through k^, and it is not necessary to evaluate B itself. 

As an example, illustrating the present method, we might consider the 
symmetric duct 


b'2(x) 


1 + Y - Y 


V/(l + 


|xir 


(13) 


for which 



where - k^. For « 1, we find 



(To O(i^) this gives s = y.) Equation (14) may be substituted into 
equations (8) and (8a) to obtain R and T for the duct (eq. ((13)). This 
example is considered further in the appendix. 

DISCUSSION 

It is noted above that, for sytnmetric ducts, when kg = 0, T « 1/2. 

The zeros of k occur at x^ and xk; thus kg = 0 corresponds to 
Xg = -X5 = 0. In this case equation (6) is equivalent to b(0) = 8(0). 

Now, when T < 1/2, x^ is real and nonzero and equation (6) determines 
B. However, when T > 1/2, no real value of Xg exists and equation (6) 
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is no longer valid. Of course, in this case there is no region in which the 
wave is cutoff; so we are considering case (1) rather than case (3) of the 
first paragraph. However, partial reflection still occurs. To continue to 
use the circular-cosh comparison equation when there is no cutoff, a crite- 
rion for determining a is required. The observation that the turning 
points, X 3 and X 5 , coalesce at a particular value of T (T = 1/2 for 
symmetric ducts), for which b and B become equal at the throat of the 
duct, suggests the appropriate criterion for determining a after the turn- 
ing points have become imaginary is that B equal b at the duct throat. 


APPENDIX - SYMMETRIC DUCTS WITH MEAN FLOW 


When the duct is syinmetric ^ the analysis of the text can be readily 
extended to include a one-dimensional gas dynamics flow through the duct. 
The density p(ex), sound speed c(ex), velocity U(ex), and Mach number 
M(ex) = U(cx)/c(€x) of that flow are expressed in terms of the duct cross- 
sectional area by well-known isentropic flow relations (sections 80 and 90 
of ref. 4). 

The propagation equation (cf. eq. (1)) is now 



(Al) 


where c_ is the upstream sound speed (so c_k is the frequency of the 
acoustic wave). In analogy to equation (3) we seek an approximate solution 
to equation (Al) of the form 


= A(cx)exp 


» /» 


M 


(1 - M^) 


dx 


XMCs(ex)/cJ (A2) 


and, analogously to equation (2), xM = -"^xm/( 1- ** ^i) with 

2 r .,2 




Substituting equation (A2) into equation (Al) gives 


4U)(S')^/(1 - M^) = x2(ex)/Cl - M^(ex)] 


(A3) 


to 0(e0). Here 




c^k^ 


C (ex)[l - M (ex)] b^(ex) 


Integrating equation (A3) gives 


s (0) 




ds 




(X) 


- M^(x) 


dx 


(A4) 


8 


in place of equation (6). Note that must be replaced by k2/l - m 2 
in expression (7) for s^ and Sjj, and x^ and x^ are now the roots of 
x^(ex) = 0. To present order of approximation, equations (8) for R and T 
remain valid (only for symmetric ducts) if k2 = k? = k^(o.)/( 1 - M£) 
and kjQ = ik^(0)i/(1 - ft). 

Equation (A4) shows that the effect of mean flow on reflection and 
transmission is to cause a Lorentz contraction of streamwise wavelengths. 

A Doppler shift also occurs because of the exponential in equation (A2), 
but this does not affect power transmission. Aside from the Doppler shift, 
tnere is no difference between upstream and downstream propagation. 

In the cases of present interest, k^ - a^/b2 is small, so the 
effect of Mach number enters primarily through the factor k^/(l - m 2) 
in icM and x^. Because M increases toward unity within the con- 
traction, X(v] (which is a function of M(x)) is more affected by flow than 
is ic|vi (which is a function of M_). In brief, the Mach number effect 
will be to decrease « in equation (12) and hence increase transmission 
through the duct. Physically, this occurs because short waves are more 
easily transmitted than are long waves; thus the Lorentz contraction of 
wavelength enables acoustic waves to pass through the duct throat. 

This phenomenon is illustrated in figure 2, where the transmission 
coefficient for the duct (eq. (13)) has been plotted as a function of Mach 
number M_. For this calculation k » 2.0, a = 1.84118 (corresponding to 
the (1,1) mode), and c = 0.3. When - 0.15, there is no cutoff 
region, so T is greater than 1/2 at all Mach numbers: T increases with 
Mach number from its value of 0.85 in the absence of flow. At the other 
values of shown in figure 1, T is zero in the absence of flow. Mean 
flow through the duct causes it to become acoustically transparent as M_ 
tends toward unity. 

The analyses of this appendix is restricted to symmetric ducts because 
in lej^i a constant Mach number of M_, has been used. Extension of this 
method to asymmetric ducts is not trivial, for a constant Mach number then 
cannot be used in the comparison equation. 
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Figure 1. - Defining sketch for acoustic transmission through a 
converging-diverging duct 



Mach number. M. 


Figure Z - Transmission coefficient as a function of Mach number for 
asvmmetric ducts with various throat radii. Throat radius, 

(l + y2)"l/2. wavenumber k-2.0; duct-mode eigenvalue a- 1. 84118; 
maximum wall slope c ■ 0. 3. 
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